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Abstract. - The traditional nuclear shell model approach is extended to include many-body forces. 
The empirical Hamiltonian with a three-body force is constructed for the identical nucleons on the 
O/7/2 shell. Manifestations of the three-body force in spectra, binding energies, seniority mixing, 
particle-hole symmetry, electromagnetic and particle transition rates are investigated. It is shown 
that in addition to the usual expansion of the valence space within the tranditional two-body shell 
model, the three-body component in the Hamiltonian can be an important part improving the 
3 ■ quality of the theoretical approach. 
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' The many-body problem is central for modern physics. A path from the understanding 

{NJ ■ of interactions between fundamental constituents to that of the diverse physics of the whole 

system is non-trivial and involves various entangled routes. Among numerous issues, the 
. questions of forces, effective or bare, their hierarchy and renormalizations are of particular 

importance. Nuclear physics is a valuable natural arena to explore this. 

Microscopic treatments provide a remarkably accurate description of light nuclei based 
on the observed, bare, two-nucleon interactions [1]. The same calculations indicate that the 
role of three-body forces is increasingly large and non-perturbative for heavier systems. The 
many-body approaches with roots in effective interactions from mean-field (one-body) to 
shell model (one and two-body) indicate a need in empirical many-body force [2-4]. ft is 
established that, regardless of ab-initio interactions, work in restricted space always gives 
rise to many-body forces; moreover, renormalizations may successfully set the hierarchy of 
importance [5,6]. 

The goal of this work is to examine three-body forces within the nuclear shell model (SM) 
approach. This includes determination of effective interaction parameters, study of hierarchy 
in strength from single-particle (s.p.) to two-body to three-body and beyond, manifestations 
in energy spectra and transitions rates, comparison with different traditional SM calculations 
and overall assessment for the need of beyond- two-body SM. Previous works in this direction 
have shown improved description of energy spectra [2, 3, 7, 8] and the significance of three- 
body monopole renormalizations [4]. 

The effective interaction Hamiltonian is a sum H/- = Yln=i where the rotationally- 
invariant n-body part is 

ff w = EE^M) J2 T LM^)T { L tm, (i) 

a/3 L M=-L 
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the isospin label is omitted here for the sake of simplicity. The operators Tj^- (a) are normal- 
ized, i.e. (0|t£?^/ {a') (a)|0) = $aa'&LL'$MM' ■> n-body creation operators coupled to a 

total angular momentum L and magnetic projection M, (a) = Cf 2 M n {ot) a\a\ . . . a\ 

where 1 is the s.p. index. The traditional SM approach is limited by the Hamiltonian 
H2 = + which is a sum of the s.p. (n = 1) and the two-body (n = 2) terms. 
In the two-body part the coefficients C\ 2 M are proportional to the Clebsch-Gordan coef- 
ficients and the index a is uniquely identified by the s.p. levels involved. For n > 2 
the index a must include additional information about the coupling scheme, the choice of 
which, in general, is not unique. In numerical work it is convenient to define the coefficients 
C^2.. [ n ( a ) an d the normalized n-body operators T^(a) using a set of orthogonal eigenstates 

\LMa) = (a)|0) of some n-particle system, see also [9]. 

Here we study identical nucleons in a single-j O/7/2 shell; see Refs. [7,10-13] for past 
works in this mass region. The best experimentally explored systems are N = 28 isotones 
starting from 48 Ca considered as a core with protons filling the O/7/2 shell and the Z=20 
40-48Q a i so topes with valence neutrons. The experimentally known states identified with 
the / 7 / 2 valence space are listed in Tab. 1. 

The three-body interactions influence nuclear masses and result in important monopole 
terms [4]. The violation of particle-hole symmetry is related to this. Within the traditional 
SM in single j-shell, with occupancy = 2j + 1, this symmetry makes the spectra of the N 
and A^-particle systems (N = Cl — N) identical, apart from a constant shift in energy. Indeed, 
the particle-hole conjugation C defined with a - m = Ca}^ m C^ 1 = (— l) J_m a 7 _ rn transforms an 
arbitrary n-body interaction into itself plus some Hamiltonian of a lower interaction-rank 
H' n _ 1 , as follows: = (—l) n H^ + H' n _ 1 . The n = 1 case corresponds to the particle 

number N = —N + fi. For the n = 2 we obtain a monopole shift 

jj(2) = H {2) + 2N ) M , M = i J2( 2L + l ) V L ] - ( 2 ) 

Any if W is proportional to N and is thus a constant of motion, which explains the particle- 
hole symmetry for the two-body Hamiltonian. The n > 3 interactions violate this symmetry 
leading to different excitation spectra of N- and A^-particle systems. The deviations from 
exact particle-hole symmetry are seen in the experimental data, Tab. 1. The excitation 
energies of v = 2 states in N — 2 system are systematically higher than those in the 
6-particle case pointing on a reduced ground state binding. 

The j = 7/2 is the largest single-j shell for which the seniority, the number v of unpaired 
nucleons, is an integral of motion for any one- and two-body interaction [10, 14]. It is estab- 
lished experimentally that seniorities are mixed [15,16]. Configurations beyond the singlc-j 
shell [8, 17-19] have been suggested to explain the effects; however, the possible presence of 
the three-body force must be addressed. In a single-j shell the pair operators Tq , T 00 and 
the particle number N form an SU(2) rotational quasispin group. The quantum numbers 
v and N are associated with this group. The invariance under quasispin rotations relates 
states of the same v but different particle numbers N. For example, excitation energies of 
v = 2 states are identical in all even-particle systems. The classification of operators ac- 
cording to quasispin leads to selection rules. The s.p. operators associated with the particle 
transfer permit seniority change Av = 1. The reactions 51 V( 3 He, d) 52 Cr and 43 Ca(d,p) 44 Ca 
show seniority mixing as the v = 4 final states are populated [16, 18]. The one-body multi- 
pole operators are quasispin scalars for odd angular momentum, and quasispin vectors for 
even. Thus, the Ml electromagnetic transitions do not change quasispin. In the mid-shell 
the quasi-vector E2 transitions between states of the same seniority are forbidden. The 
seniority mixing between the v = 2 and v = 4 pairs of 2 + and 4 + states is expected in the 
mid-shell nuclei 52 Cr and Ca. Seniority can be used to classify the many-body operators 
T^2[ and the interaction parameters. The three-body interactions mix seniorities with the 
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Table 1: States in fa/2 valence space with spin and seniority listed in the first and second columns. A 
* denotes seniority mixed states in 3B/7/2 and the seniority shown reflects the seniority assignment 
in the 2Bfa/ 2 model. In order to distinguish among the mixed states of the same spin- label (first 
column) we include an additional subscript reflecting the order in which they appear in the energy 
spectrum. Following are columns with data for iV = 28 isotones and Z = 20 isotopes. Three 
columns for each type of valence particles list name and excitation energy, experimental binding 
energy, and energy from the three-body SM calculation discussed in the text. All energies are in 
MeV's. 



exception of interaction between v = 1 nuclcon triplets given by the strength 1 . 

To determine the interaction parameters of H3 we conduct a full least-square fit to the 
data points in Tab. 1. This empirical method, which dates back to Refs. [10,20], is a part 
of the most successful SM techniques today [21]. Our procedure is similar to a two-body fit 
in sec. 3.2 of Ref. [22], but here the fit is nonlinear and requires iterations due to seniority 
mixing. In Tab. 2 the resulting parameters are listed for the proton (fixed N = 28) system 
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Table 2: Interaction parameters of 2B/ 7 / 2 and 3B/ 7 / 2 SM Hamiltonians determined with the least- 
square fit are given in keV's. 



and neutron (fixed Z = 20) system. The two columns in each case correspond to fits without 
(2B/ 7 / 2 left) and with (3B.f 7/ / 2 right) the three-body forces. The root-mean-square (RMS) 
deviation is given for each fit. The confidence limits can be inferred from variances for each 
fit parameter given within the parentheses. 

The lowering of the RMS deviation is the first evidence in support of the three-body 
forces; for Z = 28 isotones it drops from 120kcV to about 80keV. All three-body parameters 
appear to be equally important, excluding any one of them from the fit raises the RMS by 
about 10%. In contrast, inclusion of a four-body monopole force based on v = 0, L = 
operator led to no improvement. The fit parameters remain stable within quoted error-bars 
if some questionable data-points are removed. The energies resulting from the three-body fit 
are listed in Tab. 1. Based on the RMS alone this description of data is quite good, even in 
comparison with the available large-scale two-body SM calculations in the expanded model 
space [19,23]. 

The proton and neutron effective Hamiltonians are different, Tab. 2. The s.p. energies 
reflect different mean fields; and the two-body parameters, especially for higher L, highlight 
the contribution from the long range Coulomb force. However, within the error-bars the 
three-body part of the Hamiltonians appears to be the same which relates these terms to 
isospin-invariant strong force. 

A skeptic may question some experimental states included in the fit, thus we conduct 
a minimal fit considering binding energies of ground states only, similar to Ref. [10]. We 

(31 • • (3) t (2) 

include a seniority conserving part given by Vy£ with v = 1 triplet operator T jm ~ a j m ^oo ■ 
This interaction is the main three-body contribution to binding and is equivalent to a 
density-dependent pairing force [24]. In a single- j model it can be treated exactly with 

a renormalized particle-number-dependent pairing strength = + ^^E§Vj 3 ^ ■ From 
relations in Refs. [10, 13], the ground state energies with v = or 1 are 

E = eN+^^^(n-N-u)^ + (N-2 + u)M' S j, (3) 

where M' = M + j^V-^. With a linear least-square fit and Eq. (3) we determine s.p. 



p-4 



Manifestation of three-body forces in / 7 / 2 -shell nuclei. 





N = 28 


Z = 20 




-9703(40) 


-9692(40) 


-8423(51) 


-8403(55) 




-2354(80) 


-2409(110) 


-3006(120) 


-3105(156) 


M 


1196(40) 


1166(50) 


-823(55) 


-876(76) 


yd) 




18(20) 




31(31) 


RMS 


50 


46 


73 


65 



Table 3: Interaction parameters for the minimal / 7 / 2 SM determined with the linear least-squared 
fit of 8 binding energies. The variances for each parameter are shown within parentheses. The two 
columns for isotopes and isotones are fits without and with the three-body term. 



energy e, pairing Vq , monopole M, and 3-body interaction Vy J (3) using 8 binding energies. 
The results, shown in Tab 3, are consistent with the full fit in Tab. 2, the repulsive nature 
of the monopole Ty J is in agreement with other works [25] . 

Next we concentrate on 52 Cr, Fig. 1. In addition to the 2B/7/2 and 3B/ 7 / 2 interactions 
from Tab. 2 we perform a large-scale SM calculations 2B /7/2P (which includes pi/2 and 
P3/2) and 2B/p (entire /p-shell, truncated to 10 7 projected m-scheme states) using FPBP 
two-body SM Hamiltonian [26]. The 2B/7/2P model and its results are very close to the 
more restricted SM calculations in Ref. [23]. 

The seniority mixing between the neighboring 4^ and 4^ states leads to level repulsion. 
The observed energy difference of 400 keV is not reproduced by 2B/7/2 (84 keV). The dis- 
crepancy remains in some of the extended two-body models: 2B/7/2P and [23] (200 keV). 
The full 2B/p replicates the splitting, but possibly at the expense of excessive intruder 
admixtures which distort the spectrum, see also Ref. [27]. The 3B/ 7 / 2 is the best in repro- 
ducing the experimental spectrum, Fig. 1. The 3B/ 7 / 2 model predicts seniority mixing; the 
u(4f) = 2.82 and ^(4^ ) = 2.71 are inferred from the expectation value of the pair operator 

(T^Tffl) = (N - v)(2j + 3 - N - v)/(4j + 2). The 2f , however, is relatively pure with 
v{2+) = 2.006. 

Violations of the quasispin selection rules are observed in nuclei [12, 15, 16, 28, 29]. The 
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Fig. 1: Spectrum of 52 Cr. 
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Table 4: B(_E2) transition summary on 52 Cr expressed in units e 2 fm 4 , the last row being the 
quadrupole moment of the first 2 + state in units efm 2 . The data is taken from [31]. The zeros in 
the second column for the 2B/ 7 / 2 model are the results of the mid-shell seniority selection rules. In 
this model the states 2i, 2 2 , 4i, and 4 2 have seniorities 2, 4, 4, and 2, respectively, see Tab. 1. '*^In 
2E> /7/2P and 2B/p models we use 0.5 (neutron) and 1.5 (proton) effective charges. The overall 
radial scaling is fixed by B(E2,2i -> 0). (1) The life-time error-bars are used. (2) There are conflicting 
experimental results on the life-time; here we use DSAM (HI, xn7) data from Ref. [31], which is 
consistent with [30]. ^ Other independent large-scale shell model calculations obtain similar values 
for the quadrupole moment for the 2+ state: -12.3 [32], -15.0, -16.2 and -17.5 [33], all in units efm 2 . 
^ The experimental data is from [34]. 

two-body shell models beyond single- j f 7 / 2 are used to explain this seniority mixing [11, 
12,17,18,23], however there are inconsistencies, suggesting a more in-depth consideration. 
Large variations of effective charges are needed to explain electromagnetic transitions [30]. 
The particle transfer spectroscopic factors show excessive amount of components outside 
the valence space [15]. Furthermore, while going beyond the single j-shell is imperative 
to account for the observed g factors [22], the agreement between the experiment and the 
common large-scale shell model results is not perfect, Ref. [27], which may be a hint for the 
presence of a three-body force. 

In Tab. 4 B(_E2) transitions rates from the four models are compared to experiment. The 
combination of nuclear radial overlap and effective charge is normalized using the observed 
E2 rate for the transition 2i — > in the 2B/ 7 / 2 , 2B/ 7 / 2 p, and 2B/p models. The normaliza- 
tion parameter for the 3B/ 7 / 2 model is identical to the one used in 2B/ 7 / 2 - The insignificant 
difference in 2i — > B(E2) between the 3B/ 7 / 2 and 2B/ 7 / 2 shows a small admixture of v = 4 
in the 2f state. The v — 4 and v — 2 mixing in 4^~ and 4^ effects the transitions involving 
these states. For example, E2 transitions 4 2 — + 2i and 6 — * 4 2 are no longer forbidden. 

The last row in Tab. 4 compares the quadrupole moment of the first 2 + state in 52 Cr 
between models and experiment. It is quoted in units of e fm 2 . Because of the seniority se- 
lection rule the quadrupole moment is exactly zero in the 2B/ 7 / 2 model; the small seniority 
mixing between 2 + states in the 3B/ 7 / 2 model leads to the quadrupole moment of a correct 
sign but it is lower than the observed one. The results from extended two-body models in 
Tab. 4 and in the literature [32,33,35] improve the picture. However, strong configuration 
mixing is needed to generate seniority mixing which is inconsistent with g factor measure- 
ments [27] and results in a quadrupole moment higher than observed. This again invites the 
possibility of an additional three-body Hamiltonian. 

In Tab. 5 proton removal spectroscopic factors are compared between theoretical mod- 
els and experiment 51 V( 3 He,d) 52 Cr [15]. The 3B/ 7 / 2 model is good in its description of 
observation especially for the 4^ and 4^ states. It was pointed out in Ref. [15] that the 
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Table 5: Proton removal spectroscopic factors. The experimental data is taken from 51 V( 3 He,d) 52 Cr 
reaction [15]. Within error-bars of about 0.1 this data is consistent with other results [31]. 

spectroscopic factors for the 4 + states probe the v — 2 component, and, thus, their sum 
within the / 7 / 2 valence space is 4/3; this result is consistent with the observation [15] but 
does not support too much configuration mixing from outside the /7/ 2 shell which would 
reduce the spectroscopic factors and their sum. 

To conclude, the study of nuclei in 0/ 7 / 2 shell shows evidence for three-body forces. 
We extend the traditional shell model approach by including three-body forces into con- 
sideration, we find that a successful set of interaction parameters can be determined with 
an empirical fitting procedure. With a few new parameters a sizable improvement in the 
description of experimental data is obtained. The apparent hierarchy of contributions from 
one-body mean-field, to two-body, to three-body and beyond is significant; it assures the 
possibility of high precision configuration-interaction methods in restricted space and sup- 
ports ideas about renormalization of interactions. The three-body forces observed in this 
study appear to be isospin invariant. The new Hamiltonian with three-body forces describes 
well the observables that are sensitive to such forces and to seniority, and particlc-holc sym- 
metries that are violated in their presence. A good experimental agreement is obtained 
for the features of spectra, for electromagnetic E2 transition rates, and for spectroscopic 
factors. While good agreement with experiment can also be obtained by configuration mix- 
ing within some advanced large-scale two-body shell model calculations, it appears that 
the three-body effective interaction has a somewhat different manifestation which can be 
experimentally identified. In particular, with the three-body force, the seniority forbidden 
transitions become allowed without generating excessively large quadrupole moment or, in 
the case of spectroscopic factors, without loosing too much strength away from a single- j 
shell. Unlike the fit to experimental data discussed in this work, it appears to be difficult to 
fit the large-scale shell model results with a single- j Hamiltonian containing a three-body 
force. This again suggests that configuration mixing in not necessarily equivalent to a three- 
body force. The work in this direction is to be continued. It is important to conduct similar 
phcnomcnological investigations for other mass regions and model spaces; on the other side, 
renormalization techniques that would link fundamental and phenomenological forces [5] 
have to be searched for. 
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